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Abstract In a previous paper ( [physics/0203061 ) " Floating bodies of equilib- 



rium I" I have shown that there exist two-dimensional non-circular cross-sections 
of bodies of homogeous density p 7^ 1/2 which can float in any orientation in wa- 
ter, which have a p-fold rotation axis. For given p they exist for p — 2 different 
densities. However, this was found only in a Taylor-expansion in a parame- 
ter which described the distortion from the circular shape up to seventh order. 
Here a differential equation for the boundary curve is given and the non-circular 
boundary curve can be expressed in terms of an elliptic integral. 

1 Introduction 

Stanislaw Ulam asks in the Scottish Book Q (problem 19), whether a sphere is 
the only solid of uniform density which will float in water in any position. In a 
recent paperQ I considered the two-dimensional problem for p 7^ 1/2. (The case 
p = 1/2 has been solved by Auerbach ||). I was able to obtain non-circular 
two-dimensional cross-sections of bodies which can flow in any orientation in 
water by a perturbative expansion around the circular solution. These cross- 
sections have a p-fold rotation symmetry. In polar coordinates r and ip the 
boundary-curve could be expanded in powers of a deformation-parameter e 

00 

r(ip) = r (1 + 2ecos(pip) + 2 c„(e) cos(np'0)). (1) 

n=2 

The coefficients c„ = 0(e n ) were determined up to order e 7 . Although one has 
solutions for p— 2 different densities p, it turned out that r{ip) was the same for 
all these densities. 

Here a non-perturbative solution is given. It is shown that the boundary- 
curve obeys the differential equation 

. = = = ar 2 +b + cr- 2 , (2) 
yr + r' 2 
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with r' = dr/d'0- This equation can be integrated 



1p = ± 



(ar 4 + br 2 + c) 2 

dr(ar 4 + or 2 + c) 



'■\/r 2 — (ar 4 + br 2 + c) 2 

± f' 2 dg(ag 2 + feg + c) 
J 2q y / q- (aq 2 +bq + c) 2 

and is thus given by an elliptic integral. With increasing ip the radius r oscillates 
periodically between the largest and the smallest radii r> and r<, resp. The 
three constants a, b, and c are determined by these extreme radii, and by the 
periodicity of the boundary. Since r' vanishes for the extrema of r, one has 

r >,<= ar >,< + br >,< + c - ( 5 ) 
The periodicity is given by 

dr(ar 4 + br 2 + c) tt 



< T\J r 2 — (ar 4 + br 2 + c) 2 V 



(6) 



The solution of the differential equation (0) expanded in powers of e agrees 
completely with the expansion obtained in H] . 

I explain now, how I arrived at this differential equation. Denote the two 
intersections of the water-line with the cross-section by L and R. From the 
general arguments given in [^| one knows, that the length 21 of the water- 
line has to be independent of the orientation. (Here actually we will fix the 
orientation of the cross-section and rotate the water-line including the direction 
of the gravitational force.) The midpoint M of the line moves in the direction 
of the chord LR as the water-line is rotated. Suppose the coordinates of the 
intersections are given by L = (xm — 1,Um) and R = (xm + 1,Um)- If the 
line moves by an infinitesimal amount L' = (xm — I + dxL,UM + dy^), R 1 = 
(ijii + I + dxR,yM + dyji), then the condition, that the distance L'R! = 21 is 
fixed yields dx^ = dxR and the condition, that M moves in the direction of LR 
yields dyL +dya = 0. From both conditions we find, that under this change the 
same length of the perimeter du = \J dx\ + dy 2 L = dx 2 R + dy 2 R disappears on 
one side below the water-line and appears on the other side above the water-line. 
Then dxL = dx^ = cos/3dw, dyL = —dyn — sin/3du. Thus the angles between 
the tangents on the boundary and the water-line are the same on both ends. 




Figure 1: Boundary curve, water-line LR, tangents, and the angles (5. 

In the following we will use these properties in two forms: (i) as line- 
condition, requiring, that the envelope of the water-lines is the loci of the mid- 
points M, so that we obtain L and R by drawing tangents on the envelope and 
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going a constant distance I along the tangent in both directions to obtain L and 
R. (ii) as angle-condition, requiring, that moving L and R the same piece u 
along the perimeter to L' and R' , the angles @l and /3r between the tangents 
and the chord L'R' always obey j3i = (3r. Both conditions are equivalent. From 
(i) follows (ii) and from (ii) follows (i). 

In a next step consider the limit of large p. Then r will oscillate on a distance 
of order f/p, which is small in comparison to f itself. I assume, that then one can 
neglect the overall average curvature 1/f and one may find a periodic function 
/ of £ = ftp, which obeys the line-condition. Indeed assuming, that this holds 
for I 3> f/p, that is on a distance of many periods, one can derive in this linear 
case a differential equation, which has to be obeyed by the function /. This will 



be done in subsection 2.1. Expanding the solution in powers of e is in agreement 
with the terms leading in 1 /p which were found in eq. 86 of S . 

Unfortunately this procedure is not applicable for finite p. Here I work with 
a new conjecture: Although for physical reasons the periodicity p has to be an 
integer, one may try it with non-integer p. If p differs from an integer infmitesi- 
mally, then going around by the angle 2tt one arrives at the same function shifted 
however, in ^-direction by an infinitesimal angle 2a. I conjecture that there is 
also a line between L and R belonging to angles ipL,R differing by Sip w 2tt of 
constant small length 21, which obeys the line-condition. This condition yields 
a differential equation. For the linear case this is done in subsection 2.2. It is 
in agreement with that derived for large distance I. For the circular case one 
obtains the differential equation in subsectio n |3.l| , which is identical to eq. (|J). 

Finally it is shown in subsections 2.3 and 3.2, that the angle-condition holds 
for the curves described by the so obtained differential equations which are 
shifted by arbitrary distances in the linear case and by arbitrary angles in the 
circular case, resp. This last property can now also be used for two points L 
and R on the same curve, and thus shows that eq. (||) yields the boundaries for 
cross-sections with p ^ 1/2. 

The derivation for the linear approximation is given, since it is somewhat 
simpler, but one may skip section |^ and go directly to the derivation in section 

I 



2 Linear Approximation 

Consider now the above mentioned linear case. Let us start from the loci of the 
middle points M of the chord LR, described by Cartesian coordinates (<?,£)■ 
Then drawing the tangent on this curve and moving along the tangent by the 
distance I in both directions one obtains the points R and L, which should lie 
on the same curve described by /(£). Thus the functions g and / should obey 

g(fl± / g,(C) = = f(£± 1 ). ( 7 ) 
2.1 Large Distance 

Assume now, that the functions are periodic with periodicity 2 A and /(£+A) = 
— /(£) = /(—£), similarly for g. Then f(nX) = for integer n and / has extrema 
at £ = oA/2 with odd o. If we seek solutions for large I, then we put I = l a + S 
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with l a = oA/2 and small S. We write 



. " =lo + 8 u ft = , lo + d =-l = S-l g' 2 /2 + ... (8) 

Vi + fl' 2 (0 ^1 + ^(0 

and obtain by expansion in Si and addition and subtraction of the equations 
9(0 = (-) ( °- 1)/2 (<5 1 /'(e + A/2) + |/'"(e + A/2) + ...)(9) 
(lo + S) J}®. = H (o - 1)/2 (/(£ + A/2) + f/"(£ + A/2) + ...). (10) 

Note that /(£ + l a + Si) = (-) (o_1)/2 /K + I + ^i) and /(£ ~ k - &) = 
(— )(° +1 )/ 2 /(^ + -| — <5i). If / is of order i", then we realize from the second 
equation, that g' and thus g is of order 1/Z Q . The first equation yields, that Si 
is of order \jl . We put S — c/l a . Furthermore we write (— )^ 0_1 ^ 2 Z 5 = <7 and 
collect the leading terms of these two equations 

g=(c-g' 2 /2)f, g' = f, (11) 

where the argument of g is £ and of / is £ + A/2. From these two equations we 
obtain differential equations for / and g. Taking the derivative of the second 
equation and substituting into the first one yields 

g=(c-g' 2 /2)g". (12) 

If we multiply this equation by g', we can integrate it and obtain 

g 2 = (c-g' 2 /4)g' 2 + c 2 . (13) 

In order to obtain a differential equation for / we solve this equation for g and 
differentiate it. This yields 

t - fMz£m ,i4) 

and substituting g' = f the equation 

Vc2 + cf 2 -fy4=(c-f 2 /2)f (15) 

or 

(c-/ 2 /2) 2 (l + /' 2 )=c 2 + c 2 , (16) 
which we will use in the form 

; 1 - af + b (17) 
with new integration constants a and b. 
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2.2 Infinitesimal distance 

We perform now a similar calculation, but assume, that we can represent /(£+ct) 
and /(£ — a) similarly by g, where now a and I are infinitesimally small. Thus 
we assume 

g(Q± /^ = = /((Tai ■ 1 )■ (18) 
Vi + s' 2 (0 Vl + 5 /2 (0 

It will turn out, that we need the expansion in a to first order and in / in third 
order. We obtain 

9 = f+ oM I'm /'' (19) 



2(1 + g' 2 )- 



l 2 



9 = / / - T / / V / T+F+ g/1 , (20) 



Z J v 3 6(1 + 5' 
Differentiating the first equation yields 



,_ fl , i 2 f" i 2 9'g"f" r91 s 

5 + 2(l + . 9 '2) (1 + ^)2- ^ 

Equating both expressions and choosing a = cl 3 yields 
, fin „i„n fii 

c/ '^^ + 5(iWtiT7? = - (22) 

In the limit a — > one has g — f ■ With this substitution and multiplication by 
1 / \/l + /' 2 we ma y integrate this equation 

c/+ 3(l + fy +Cl = °' ^ 
Multiplication by /' allows a second integration which yields 

7r = = |/ 2 +3c 1 / + 3c, (24) 

Comparison with the equation obtained for large distances ( |l"7| ) shows agree- 
ment, if we change the notation for the constants. There is no linear term in 
eg-dTif) like the term 3ci/ in eq.(24). The reason is, that in eq. ( |l7| ) we required 
/ to vanish on the average. But of course a constant can be added to / and g, 
yielding such a linear term. Thus both procedures yield the same differential 
equation. 



2.3 Arbitrary Distance 

We show now the following theorem. Consider two curves Cl and Cr which 
are governed by eq. (|l^) with the same constants a and &, and point L € Cl 
with Cartesian coordinates £l, Jl and R G Cr with coordinates £r, Jr. If now 
the angle /3l between the tangent in L and LR equals the angle (3r between the 
tangent in R and RL, and Jl fa, then this property will hold between any 
two points L 1 e Cl and R' 6 Cr, which are obtained by moving from L and R 
by the same arc u in the same direction along the curves. Simultaneously the 
distance L'R! = LR is constant. 
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In order to show this we first express the angles (3l and f3 R . Let us denote 
the angles of the curves against the £-axis by <Pl,r, 

tan^fl = f' LR , cos<j) LM = 1 =afl R + b, (25) 



sm<j> LM = = ^/l - {gfl + 6)2 =; Wl ,r (26) 

V 1 + /£* V 



and the angle of LR against the £-axis by 7, 



tan 7 =|-|, - 7 =^A, 0067 = ^1-^=^=:^, (27) 
where 21 is the distance Li?. 




Figure 2: In this figure 4>l is given by /?l + 7. The angle is negative in this 
figure, and (3 R = 7 - tj> R . 



Then we have Pl — 4>l — 7, Pr = 7 — <Ar- The condition /?£ = /3r can now 
be expressed in various ways. From cos(</»l — 7) = cos(7 — <f> R ) one obtains 



(a/| + b)W. t + ^-/^W L = (af 2 R + b)W 7 + f -^J±W h 



"> ' 21 
which for f R ^ Jl reduces to 



21 



We show now, that 



a(f R + f L )W 7 



n:= (3 L ~ F 



W L -W R 
21 



^-2 7 



(28) 



(29) 



(30) 



is independent of u, if initially it vanishes, n = 0, and f R ^ /l holds. For this 
purpose we calculate the derivative ^ and show, that it vanishes, and thus, 
that n remains 0. For this purpose we calculate the derivatives needed and 
indicate the derivative with respect to u by a dot, 



Il,r — 
W L , R = 



sin((/)i i fl) = W l ,r, 

-T? Jl,r = -2af LtR {af LR + b), 

ajL,R 

(fL-f R )(W L -W R ) 



(31) 
(32) 

(33) 
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We use generally for the derivatives of the angles 



= cos (j)(sm (j>) — sin </>(cos (f>) (34) 

and obtain 

4>lm. = -2af L ,R (35) 

and 

1 = W ^— Ww, • (36) 

We make now use of (|29| ) and obtain 

j = -a(h + fa), (37) 

from which we conclude h = 0, provided /l ^ f R ■ 

For Jl = Jr one has Wl = ±Wr. If the minus-sign applies, then n vanishes 
for all Jl and I. In this case, however, W 1 — 1, W 7 — and we obtain 
h = —4af + 2Wl/1- One can easily check, that for f R = /z + e one obtains 
n = e(2af /Wl — l/l) + 0(e 2 ), so that n — holds also as we approach Jl = Jr- 
Thus the theorem is proven. 

Obviously instead of considering two curves, we may consider two points L 
and R on the same curve, which implies, that the midpoints of these chords are 
the envelope of these lines as we increase u. 



3 Circular Case 

The point R with polar coordinates (r, ip + dip) lie on the tangent of the envelope 
through M with polar coordinates (/Lt, a distance / apart. Then we have 

r 2 =fl 2 + -Wl + { 2 (38) 

Denote the angle OMR by 7. Then 7 is related to the slope of the tangent by 
— cos 7 = n' / \J fj? + fi' 2 and sin 7 = /1/ -\//i 2 + [i! 2 ■ For Sip we obtain 
tan Sip = lsm~f/(fi — I cos 7) = Z/i/(/xi/ fi 2 + fi' 2 + Z/i'). Thus we obtain 



.//i 2 ±^^= + P = r(^Ta±arctan ^ ), (39) 

where the upper sign yields on the curve rotated by a and the lower sign tl 
on the curve rotated by —a. 

3.1 Infinitesimal angle 

We expand both sides of eq. ( |39"| ) up to order I 3 and a and add and subtract 
the equations for both signs. Then we obtain 

At+ 2(^ 2 + M ' 2 ) r+ ( ~^ 2 + M ' 2 ) _ 2(/i 2 + /i' 2 ) j ' 1 j 

*y _ , q / / 2 , /2 1 , 2 /(-m 2 +3^ 2 ) 

2(/i 2 + M ' 2 ) r ; r VM+Ai 3^ 2 (/Li 2 + /z' 2 ) 

r'V r'" 

:)• (41) 



fi(fi 2 +H 12 ) 6(^ 2 + M ' 2 )' 



7 




Figure 3: The angles in the circular case. The angles <f> determined by tan <fi = 
r'/r are related to (3 and 7 by <Pl — 7l + Pl — 7r /2, ~4>r — Jr + (3r — n/2. 

Now we differentiate eq. ( |40| ) 

/ , ,2/m'(-m 2 + m /2) mmV n 



2( M 2 + /i' 2 ) 2 (^ 2 + /1' 2 ) 2 y 
, a rV 2 (3/x 2 + /x /2 ) _ r'V (2// 2 + 

1 M 2 Gu 2 +// 2 ) 2 fi(n 2 +ii' 2 ) 2 
fj,"r'(-fj, 2 + fi' 2 ) r»n"it' r" 1 



' /i( M 2 +M' 2 ) 2 (M 2 +M' 2 ) 2 ' 2( M 2 + M ' 2 ) ) (42) 

and subtract eq. (|4l| ) from this derivative. The resulting equation contains 
terms proportional to a/l and to I 2 . Thus we set a proportional In the limit 
I — > the function [i approaches r and similarly for the derivatives. Thus we 
have 



/: = --^r'v / r 2 + 
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I 3 

r'(r 2 + 4r' 2 ) r ' r "(r 2 - r' 2 ) 



3(7*^ _j_ j,^2 _|_ ^/2^2 ^j.2 ^,12^2 3(7* 2 ~|~ 7* /2 ) 

This differential equation can be integrated. One finds, that 



0.(43) 



with 



_ Vr 2 + r' 2 d/i _ Vr 2 + r' 2 df 3 
r dtp r 3 dtp 



(45) 



3(r 2 + r' 2 ) 3 / 2 2/ 3 ' 



h = r 3 (r + r") _or^ 

J3 3 ( r 2 +r /2)3/2 4 p ' V ») 

Thus /1 and f 3 are constants. Elimination of r" yields the differential equa- 
tion of the bounding curve quoted in the introduction, eq. (|[). 

' =ar 2 + b + cr~ 2 (47) 



y/r 2 + 

with constants a = \im3a/l 3 , b = 3/i/2, and c = —3/2/2 



3.2 Arbitrary angle 

Now we consider the two curves rotated against each other by an arbitrary angle 
2a and perform similar calculations to those for the linear case. 
We hrst express the angles 

0L = <t>L ~ 1L + J, 0R = ^-4>R- 1R- (48) 

and show the following theorem: If the condition 0l = Pr, that is 

n := 4>l - Jl + (I>r + 1r = (49) 

is fulfilled and tl ^ r R and the points L and R move along their curves by the 
same arc u to points L' and R' , then the condition n = remains fulfilled, that 
is for all these chords (3l = f3 R holds. Simultaneously the distance LR = L'R' 
remains constant. 

In order to show this we calculate dn/du. For this purpose we first list the 
sines and cosines of the angles 

COS<f) L ,R = L,R _ _ rL R ^ aJ ,2 + b + _ j j-gQ^ 



L,R ' ' L,R 



L,R 



1 L R 

sin0 L ,fl = ' = =: W l ,r (51) 



r l,R + r 'l,R 



cos -y L ,R = — - (52) 

4:lrL,R 

sin7L,_R = (53) 
Wl,r 



W 1 = y'-16« 4 + 8P(r| +r 2 R ) - [r\ - r|) 2 (54) 
Their derivatives with respect to u indicated by a dot are 

rL,R = sin</>L,R = W LtR (55) 
w a dcos0 L , fl 

W l ,r = -cos(j) LM — — (56) 

dTL,R 

Performing the derivatives for and r R and using ( J5q ) one obtains 

Wj = (W R cos 1l + W L cos 7B ) (57) 

V V /y 

Thus we obtain using again eq. (|34|) 

• 2 dcOS0 Ljfl 2 dcOS^L,fl 

0L,fl = -cos z (j) L}R — ■ W£ R — '— 

dr L . R ' dr L . R 

dCOS(j) L:R 2 , , C /ro\ 

cos7l._rW 7 , W^(r R . L W R . L ~ r L , R W L , R ) 
1l,r = — tt— 1 ■ • ( 59 ) 
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Now we use n = in order to simplify the derivatives 7. First we evaluate 
sm(4>L + Tfl) = sin(7£ — <f>n) which yields 

w , w W-y . COS <f)R cos 

VF i? cos7L + VK L cos7 fl = — -4 ) (60) 

4/ r\L tr 

Inserting this expression into eq.(p7[) yields 

W n = 2(r R cos cf) R - r L cos <f> L ) (61) 

Further we determine the products W 7 WL t R from cos(0l — 7l) = cos(^>r + 7_r) 
and cos(4>l + In) = cos(4>r — 7l), which yield 

cos0 l cos7lH = cos0 fl cos7_R f , (62) 

AItl Wr 

W 7 W L W 7 W R 

cos 4> L cos 7 fl -\ = cos0 fl cos7 L H -± . 63) 

AlrR AItl 

For tl ^ tr one obtains 

W,W L = 8al 2 r L r 2 R + (or| + br L )(4l 2 + r| - r|) - ^ + ^ ~ ^ ,(64) 

= -8a^r| - (ar| + br fi )(4^ 2 + r| - r|) + c(4 ^ + ^ ~ ^ 65) 

Tr 

and thus 

7L = -a(r 2 L + 2r 2 R )-b+^, (66) 
7« = a(2r 2 L +r R ) + b-^, (67) 

from which one concludes h = 4>l — jr, + <t>R + 7-R = provided r/, 7^ r/j ■ 

For = t\r one has — 1r and 4>l — ±<f>R. If the minus sign applies, 
then = and n = 0. In this case we are not allowed to divide by — r R , 
which we did in our derivation. It turns out, however, that in this limit W 1 — 
and one obtains 

^ ^ + c h), (68) 



2r L l 



2 



c (r L ) = -6ar 2 L - 2b 



2c 



„2 



2 4({rl + rl)(r% + r\) - 4.1) + 2{r > T < ^ 
r Z r£(r>+r<) 



3.3 Application to the boundary 

We may now apply the theorem to our original problem. In this problem the 
two curves are identical and constitute the boundary. Suppose we have now 
a curve obeying the differential equation (|^) with integer p. Then we have to 
look for points L, R on this curve which have either the properties n = and 
r L 7^ tr or otherwise the properties = tr and n — 0, eq. (pq). We will now 
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show, that there are at least p — 2 solutions, which obey the second condition. 
For this purpose we divide the curve into pieces A\, ... A p in which r decreases 
from r> to r < as ip increases and into pieces B%, ... B p , where r increases from 
r< to r>. We now look for solutions tl — tr and h = 0, where tr is in A\ and 
tl in Bk- Thus we have to find a solution of (pf). We realize, that with the 
exception of / = the first term in (|6^ ) vanishes for r = r> and r = r<, since 
= 0, because r' = at these extrema. 




Figure 4: Boundary for p = 7 and e = 1/49 (eq. |l]) with the pieces A n and -B„. 
The outermost curve is the boundary, the inner curves are the envelopes of the 
water-lines. The innermost one for density p = 1/2, the two others for densities 
0.22748 and 0.77252, and for 0.05269 and 0.94731. 



On the other hand we can evaluate cq for these extrema 

/ x 2(r > -r < )(l + 2ar > (r > +r < f) 

co(r>) = -, ■ r (70) 

r>(r> + r<) 

/ v 2(r> -r < )(l + 2ar < (r > +r < ) 2 ) 
c o(r<) = -, j r (71) 

r<(r> +r<) 

The denominator y/r 2 — (ar 4 + br 2 + c) 2 of the integral vanishes only at the 
extreme points r = r > and r = ?*<. If we express b and c in terms of a, r>, and 
r<, then 

.2 ^„„4 , ,2 , x2 



r 



(ar + 6r + c) = (r> — r)(r — r<) I h a(r + r>)(r + r<) 



x (r> + r)(r + r<) ( h a(r - r>)(r — r<) 

™ > + J"< y 

(72) 
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Thus l + a(r + r > )(r + r < )(r > +r < ) > 0, and therefor the factors (l + 2ar>(r> + 
r<) 2 ) and (1 + 2ar<(r> + r<) 2 ) in the expressions for Co are positive. Therefore 
one has co(r>) < and Co(r<) > 0. Therefore h changes sign between r = r > 
and r = r<. Since n is a continuous function of rj,, there must be a zero in 
between. This argument holds for all pieces Bk with the exception of those 
adjacent to A±. (For those adjacent I may become 0. In this limit one finds h = 
0.) Therefore we have at least p — 2 solutions, which is in agreement with those 
found in We leave aside the question, whether for sufficiently large ratios 
r> jr < more solutions appear and whether there are other solutions for p ^ 1/2. 
We finally mention that in the near-circular case (r> — r<)/ (r> + r< ) <C 1 the 
periodicity p is related to a by p 2 = 1 + 8ar 3 , where r is some intermediate 
radius. 

4 Conclusion 

We have given a closed representation of the boundaries of logs floating in any 
direction with p ^ 1/2 in terms of a differential equation which can be solved 
by an elliptic integral. The differential equation could be obtained, since points 
of two of these curves rotated against each other can be connected by chords so 
that by progressing on the two curves by the same length of the perimeter the 
length of the chord does not change. Using that this holds also for two curves 
rotated against each other by an infinitesimal angle connected by a chord of 
infinitesimal length the differential equation could be derived. The differential 
equation of the linear analogue of this curve has also been given. 
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